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Abstract 

A rational triangle has rational edge-lengths and area; a rational tetrahe- 
dron has rational faces and volume; either is Heronian when its edge-lengths are 
integer, and proper when its content is nonzero. 

A variant proof is given, via complex number GCD, of the previously known 
result that any Heronian triangle may be embedded in the Cartesian lattice 
Z^; it is then shown that, for a proper triangle, such an embedding is unique 
modulo lattice isometry; finally the method is extended via quaternion GCD to 
tetrahedra in Z"^, where uniqueness no longer obtains, and embeddings also exist 
which are unobtainable by this construction. 

The requisite complex and quaternionic number theoretic background is sum- 
marised beforehand. Subsequent sections engage with subsidiary implementation 
issues: initial rational embedding, canonical reduction, exhaustive search for em- 
beddings additional to those yielded via GCD; and illustrative numerical exam- 
ples are provided. 

A counter-example shows that this approach must fail in higher dimensional 
space. Finally alternative approaches by other authors are summarised. 

Keywords: triangle, tetrahedron, lattice embedding, axial pose, Cayley-Menger de- 
terminant, quaternion arithmetic, Euclidean algorithm 
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1 Introduction: Rational and Heronian Simplices 

A rational triangle, tetrahedron, or simplex in n-space R" is a polytope with n + \ 
vertices, and rational edge lengths, face areas, solid volumes, etc. When the edges have 
integer length, the simplex is Heronian; when the GCD (greatest common divisor) of 
those lengths equals unity, it is primitive (Heronian). When its n-dimensional content 
is nonzero, it is proper — we generally have only proper cases in mind, despite which 
much of the following may easily be seen to apply more widely. Obviously any rational 
simplex may be dilated to a Heronian simplex. 

Our purpose here is to show that for n = 2, 3 the simplex in M" may be posed 
first axially in Q" (which is elementary), and finally embedded in the lattice Z"; 
that is, congruently so that its vertices have Cartesian coordinates not just rational, 
but integer. For triangles the result is only about 10 years old [T3]; we propose a 
construction based on complex- number GCD, and related to the approach of [H]. 

For tetrahedra the paradigm must be extended to embrace quaternion GCD, where 
commutativity and unique factorisation are no longer available: the correspondence 
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between the two arguments is illuminating, and we emphasise the analogies and con- 
trasts by presenting the main theorems in sections IH |8]in copycat format. The level 
of detail in their proofs might be felt excessive, to the point of obscuring the wood 
for the trees: it was motivated by a desire to avoid elementary blunders which have 
compromised numerous earlier versions. 

Incidentally, an obscure but significant step in the argument involves first estab- 
lishing that the preliminary rational pose (on which GCD is to act) has denominators 
with prime factors p = 1 (mod 4) only. 

It is natural to enquire how these results might generalise to higher dimensions. 
There are currently several obstructions: the quaternion representation employed can- 
not be extended further; even if it could, the theorem in the form given here would no 
longer hold; and in any case, no examples of Heronian pentatopes are known to which 
any such theorem might be applicable. 

Despite having spent much time contemplating this topic, this author still finds the 
connection revealed, between the number theoretic concept of GCD and the geometric 
concept of lattice embedding, both surprising and mysterious. 

The contribution to this project from Warren D. Smith, in the form of copious 
literature references, enthusiastic criticism, and inexhaustible energy, is gratefully ac- 
knowledged; also some motivating ideas, particularly behind the proof of the embed- 
ding theorem, originated with Michael Reid. 



2 Edges, Triangles, Tetrahedra and their Content 

The Hero formula giving the area c? of a triangle with edge lengths u,v,w as 

(4c?)^ = {u + V + w){u + V — w){u — V + w){—u + V + w); (1) 

is case ti = 2 of the Cayley-Menger formula for the content of a simplex in M", which 
deserves to be better known [T2l: 
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For 71 = 3 the volume e of the tetrahedron with edge lengths u, . . . ^ z is given by 



2(12e)2 
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Since these coeffcients turn out to be even, (12e)^ is again a homogeneous symmetric 
polynomial with integer coefficients, but now cubic in the edge squares. 

To specify a particular free tetrahedron, a correspondence must be established 
between vertices and edges: denoting vertices P, Q, i?, 5, edge-lengths will be specified 
by a hexad of integers presented in sequential order 

[gP, PP, Pg, SP, SQ, SR] . 

The amount of data output can be reduced by restriction to primitive cases, in 
which the GCD (greatest common divisor) of the edges equals unity; also by reduction 
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to canonical form, with vertices permuted so that the vector of edge-lengths represents 
the largest possible 'number', where individual components are regarded as 'digits'. 
[A humbling exercise in algorithmic design is to sort such a hexad into canonical form 
efficiently, without first generating every vertex permutation: our best effort achieved 
a maximum 22 transpositions. And we are almost sure it actually works . . .] 

Methods for constructing Heronian simplices are discussed in detail elsewhere [9], 

El, [n]. 

Rational simplices are of interest as test-cases since many of their subsidiary el- 
ements are likewise rational — for triangles, the edges, altitudes, medians, centroid, 
orthocentre, incentre, circumcentre, etc. Less obvious is the rationality of the angular 
quantities 

cotp/2 = (s — w)s/d 

etc., where p denotes the angle at vertex P. Conway christens this property "geodetic": 
it implies that all the trigonometric functions of the angles are also rational. 

3 Representation of and 1? via Complex Num- 
bers 

With x, y real in ]R, denote a complex number in C = K[t] by 

Z — X + by, 

its norm by 

||Z||=x2+y2^ 

its conjugate by 

Z = X — ty. 

Any P — X + Ly represents a finite point with Cartesian coordinate {x,y). The 
representation is inhomogeneous: scalar multiples of P represent distinct points. The 
origin point is represented by AT = 0. The distance squared between points P, Q equals 
||P-Q||. 

Any unit U = cos t + l sin t represents rotation through angle t around the origin; 
non-unit U with \\U\\ ^ 1 do not represent rotations. The identity rotation is repre- 
sented by [/ = 1. Composition of rotations is represented by commutative complex 
product. 

The transformation 

P UP where U = cos t + l sin t 

rotates the point P = x + cy through angle t. 

A (Gaussian) integer complex number lies in its components rational integers. 
A complex prime is either 2, a rational prime p = —1 (mod 4), or of form x±Ly where 
x'^ + y'^ = p with p = 1 (mod 4). Factorisation is essentially unique modulo association 
— product with unit integers ±1, ±t. 

Complex integers are a Euclidean domain: given X,Y G their greatest com- 
mon divisor GCD(A', Y) is unique modulo association. Algorithmic details of CCD 
are discussed in section |6l For further detail see [4 chap. 2. 

Assertion 1. If Z G Z[t] is primitive, then GCD(Z, Z) = 1 (or associate). 
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Proof. 

GCD(x + Ly,x- iy) = GCD(a; + ty, 2x) = GCD(a; + by, 2/,y), 

which divides the complex prime 2 = 2 GCD(a;, y); since one of a;, y is odd, the GCD 
cannot be 2, so it must be (an associate of) unity. □ 

Unit integers represent rotations through kn/2 around the origin. The full group 
of lattice isometrics is generated by these 4 rotations, reflections in axes (complex 
conjugacy etc), and translations (integer addition). 

4 Lattice Embedding Theorem in 

Assertion 2. Suppose all hut one of a set {Q} of k + 1 points already embedded in 
; the other point R rational in with common denominator r; and the squares of 

2 

all their point-to-point distances in Z. Then rotation via complex product Q ^ X Q 
where X — GCD(i?, r) embeds the entire set congruently in 1? . 

Given the first assertion, translating some point of the set to the origin and applying 
induction on the number fc of points, immediately 

Assertion 3. // a set of points lies in Q^, and their squared distances lie in Z, then 
the entire set may be embedded congruently in 1? . 

And taking the set to be the vertices of a Heronian triangle in some rational pose 
(section [7]) , with edge- lengths themselves integer, 

Assertion 4. Any Heronian triangle may be embedded congruently in 1? . 

Proof, of assertion [21 Let the points P,Q,... lie in Z^, except for R in with 
Cartesian coordinate {x/r,y/r) where r = hGD{x,y) (least common denominator). 
Represented by complex numbers, all are integer except for R; now dilating the entire 
set by r yields also R" — rR integer. 

The first lattice point P of the set may be translated to the origin P = 0; now if 
k = 1 the result becomes trivial. Otherwise since all distances squared ||P — are 
integer, all dilated points have magnitude divisible by r^. Set X = GCD(i?",r), for 
which 

r = XX= \\X\\ =XX. 
The rational point R is transformed by X to the lattice: 

I I ||P"|| = R"W, 

and by assertionfl] i?", R" have common factor GCD(x, y), which is scalar and coprime 
to r. Therefore X^\R"] and R is transformed to 

X'i?7r2 ^ x'x^P'/r^ ^ R' 

for some integer R' = X Rjr ^Tl? . 

Also each lattice point Q is transformed by X to the lattice: 

I II rg - P" II = (rg - P")(rQ - i?") 

which is transformed to the equal 

||XVg-X^P"||/r2 = (x\q -X^ R"){X\Q - X'^W)/r^; 
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hence 

I (X^rQ -X^R"){X^rQ - X^W). 
Now for r e Z, F e Z[<,] we have r | F iff r | F; so via unique factorisation 

I X^rQ-X^R"; also | x'^R", 

and adding, finally 

I rX^Q. 

Therefore Q is transformed to Q' = X^i?/r e Z^. □ 

The construction is in fact reversible, although this seems difficult to establish 
directly. The following argument relies heavily on unique factorisation. 

Assertion 5. Lattice embedding of a proper Heronian triangle is unique modulo lattice 
isometry (association, integer addition, conjugacy). 

Proof. Suppose P, Q, R and P, S. T arc; lattice embeddings into Z[l] from vertices of 
the same free triangle, translated and reflected to share the origin P = and to have 
anticlockwise sense. 

Let Q' = Q/GCD(Q,5'), S" = S/GCB{Q,S), and set U = S/Q = S'/Q'; the 
isometry U rotates P,Q,R to P,S,T. Since ||Q|| = \\S\\ we have ||Q'|| = = q 
say. 

If g = 1 then U G Z[/,] is a lattice isometry, so the embeddings arc essentially 
equivalent. Otherwise, each rational prime factor of q splits into conjugate primes; so 
Q' = S' and t/ = S'/Q' = S''^/q, and similarly U = T'^/q. But now S' = ±T', so the 
triangle was after all improper. □ 

5 Representation of and I? via Quaternions 

With s,p, q, r real in M, 

X = s+pi + gj+rk 
denotes a quaternion in H = M[i, j, k]; 

II v-ll 2 I 2 I 2 1 2 

II A II = s +p + q +r 

denotes its norm; 

X = s — pi — gj — rk 

denotes its conjugate. 

In general X ^ represents the finite point at distances p/s,q/s,r/s from the 
Cartesian coordinate planes. The representation is homogeneous: any nonzero scalar 
multiple of X represents the same point. The origin point is represented by X = 1 
(or any nonzero scalar) . A pure quaternion with scalar component s = represents a 
point at projective infinity. The distance squared between points P, Q equals ||P — 
only when both have been co-normaliscd, with scalar parts unity. 

Dually X represents rotation through angle t around a line meeting the origin, 
where s = cotf/2 and the direction cosines (of angles made by the line with the 
coordinate planes) are p, q, r. The representation is homogeneous: any nonzero scalar 
multiple of X represents the same rotation. The identity rotation is represented by 
X = 1 {or any nonzero scalar). Composition of rotations is represented by quaternion 
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product, which is in general non-commutative, unless its arguments differ only in their 
scalar component. 

In vector terms, given 3-vector u, unit 3-vector v, and angle t, the transformation 
(group-theoretic conjugation) 

P YPX where X ^ cott/2 + v 

rotates the point P = 1 + u through angle t about axis v. 

Interpreted as rotations, quaternions are analogous to conventional polar or con- 
travariant vectors; as points, to axial or covariant vectors: and this duality accounts for 
their varying normalisations. Again, transformation by a rotation is the only situation 
in which it is legitimate to multiply a point. 

Symmetries of the Cartesian lattice are generated by translations (integer addi- 
tions), reflection in the origin (conjugacy), and 24 rotational fixing the unit cube via 
transformations: 

1 (identity ) x 1; 

i, j, k (tt around axis ) x 3; 

(l±i)/V2, ... (7r/2 around axis ) X 6; 

( i ± j)/\/2, ... (''■/2 around edge perpendicular ) x 6; 

(1 ± i ± j ± k)/2 (27r/3 around diagonal ) x 8. 

A (Lipschitzian) integer quaternion G j, k] has rational integer components. It 
is primitive when the GCD of its components equals unity. Rational points are repre- 
sentable by integer quaternions, lattice points by integers with unit scalar component; 
and conversely. 

Quaternion number theory is complicated by non-commutativity of multiplication 
and non-unique factorisation. A prime X € j, k] has prime norm ||^|| = p, while 
X ^ Z itself; an integer quaternion may factorise into primes in several essentially 
distinct ways. Factorisation is discussed in detail in |3] chap. 5. 

Integer quaternions are almost a Euclidean domain: via [5| Theorem 3, provided at 
least one of Y, Z e Z[ i, j, k] has odd norm, their left- and right-hand greatest common 
divisors 

GCDL(y,Z), GCB r{Y,Z) 

are guaranteed to exist — though in general distinct — and are unique modulo asso- 
ciation — product with units ±1, ± i, ± j, ± k. 

[However provided q \ Q and q E Z, then via [6] Theorem 2 there do still exist 
integers X, W such that q = XX and Q = XW; but only when q is odd can X = 
GCDL(Q,g) be guaranteed unique. Analogous observations apply on the right-hand 
side.] 



6 Algorithms for GCD 

Define rounding by {x) = [x + 1/2J G Z, the integer nearest to x G R; and for 
quaternions 

(X) = {s+pi + qj + rk) = (s) + (g)j + (r) k; 
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then left and right remainder functions are defined by 



XmodLY = X - Y{YX/\\Y\\), 
XmodnY = X- {XY/\\Y\\)Y; 

where — Y~^ equals the commutative multiplicative inverse of Y. 

Now the Euclidean algorithm can be employed, enhanced to catch those cases where 
no GCD exists: 

GCDi(y,Z)^ 
while Z 7^ do 

X ^Y; Y ^ Z- Z ^ X modL F; 

if ||Z|| > ||y|| then ABORT fi od; 
return Y\ 

For GCDj^(y, Z), substitute modR for modL above. 

This procedure is guaranteed to terminate only provided at least one of ||y|| or 
||z|| is odd, see [5] Theorem 3. For example when y = 2, Z = l+ i+ j+ k (both 
with norm 4), the distinct primes 1 + i, 1 + j, 1 + k all divide both, but no product 
of all three (with norm 8) can possibly divide either; so no common divisor X can be 
divisible by every divisor of both. 

GCD in C follows lines similar to GCD in H, though simplified by commutativity, 
just two components, and existence for all arguments. 

7 Rational Embedding of Heronian Tetrahedra 

Assertion 6. In order that a simplex in M" be rationally emheddable in Q", it suffices 
that the edge-lengths have rational squares, and there exists some flag (incrementally 
ascending sequence of vertex subsets) all having rational content. 

The result extends in a natural fashion to polytopes with more than n + 1 vertices, 
see |8j Prop. 3.3; here we are principally concerned only with the case n = 3. 

Proof. Given an arbitrary free tetrahedron with edge lengths u, v, w, x, y, z to be posed 
in Q^, its vertices may be posed in axial subspaces of increasing dimension thus: 

P = l, Q = l + Qii, i?= l + i?ii + i?2j, 5 = l + 5ii + 52j + 53k. 

Elimination of the 6 equations ||P — (5|| = w etc. displays the vertex coordinates as 
rational functions of edge length u of PQ, face area d of PQR, volume e of PQRS, 
and the squares of other edge lengths: 

Qi = w, 

Ri^{v^ -w^ + Qi^)/2Qi, 
i?2 = 2d/u, 

Si = {x^-y^ + Qi^)/2Qi, 

5*2 = {x^ -z^ + Rf^ + i?2^ - 2 Si Ri)/2 i?2, 

5*3 — 3e/d. 

□ 
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Each vertex will then be rescaled homogeneously as a primitive integer quaternion, 
with scalar component the LCD of the reduced rational coordinate components, for 
submission to the lattice embedding procedure in section |8l [The expressions are only 
valid for length u and area d nonzero: a more elaborate algorithm is required to cope 
with all improper cases.] 

The denominators of these rationals possess an obscure property, unexpectedly nec- 
essary to permit the use of GCD in construction of lattice embedding transformations. 

Assertion 7. Each coordinate component of each vertex of a Heronian tetrahedron in 
axial pose has denominator the product of primes p = 1 (mod 4). 

Proof. Borrowing a technique from the coordinate components are expressed as 
homogeneous rationals in terms of cotangent half-angles: at triangular vertices in 
R^, and at dihedral edges in R'^. These cotangents are rationals /'//" ^ Q with 
GCD(/', /") — 1; when substituted into the coordinates, they give rise to denominators 
with factors of form /'^ -|- /"^ . 

A theorem from elementary number theory due to Euler [Tj guarantees that f'^+f"^ 
factorises into primes p — Ak+\, possibly with a single extra factor 2 when both /', /" 
are odd. It then remains only to check that these surplus 2's cancel with the numerator 
in each component. 

So consider the vertex 5" above, which subsumes all other cases. Denote by 
k = 2c/u the height of face PSQ; at vertex P we find g in terms of S etc. 

cosP = S'i/a:, sinP = fc/.T, 
5 = cotP/2= (1 + cosP)/sinP = {x + Si)/k; 

and similarly h at dihedral U along edge PQ 

cosU = S2/k, sinC/ = S's/fc, 
/i = cot [7/2 = (1 + cost/)/ sin C/ = {2 c/u + S2)/ S3; 

then inverting to get S in terms of g, h etc. 

S = 1 + X cos Pi + X sin P cos U j + x sin P sin U k 
= l + {g'^l)/ig^ + l)-xi 

+ 2g{h'-l)/{g^ + l)ih^ + l)-xi 

+ A gh/ {g^ + l){h^ + l)-xk 
= l + ig''-g"')/ig''+g"')-xi 

+ 2g'g"{h'^ h"^)/{g'^ + g"^){h'^ + h"^) ■ 

+ Ag'g"h'h"/{g'^+g"^){h'^ + h"^)-x\., 

where g = g' / g" and h — h' /h" are positive reduced fractions. If it is even, (/'^+/"^) = 
2 (mod 4) and (/'^ — /"^) is also even; so all 2's cancel from the denominators as 
required. □ 

8 Lattice Embedding Theorem in 'E? 

Assertion 8. Suppose all hut one of a set {Q} of k + 1 points already embedded in 
1? ; the other point S rational in Q'^ with common denominator s; and the squares of 
all their point-to-point distances in Z. Then rotation via quaternion transformation 
Q — >■ XQX where X = GCD]^{S, s) embeds the entire set congruently in 1? . 
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Given the first assertion, translating some point of tlie set to tlie origin, and ap- 
plying induction on the number k of points, immediately 

Assertion 9. // a set of points lies in Q'^, and their squared distances lie in Z, then 
the entire set may be embedded congruently in 1? . 

And taking the set to be the vertices of a Heronian tetrahedron in some rational 
pose (section [7]) with edge-lengths themselves integer. 

Assertion 10. Any Heronian tetrahedron may be embedded congruently in 1? . 

Proof, of assertion [HI Let the points P,Q,... lie in Z'^, except for S in with 
Cartesian coordinate {x/s,y/s, z/s) where s = LCD(x,?/,z). Represented by integer 
quaternions, every point is primitive; and each has scalar component unity, excepting 
S = s + xi + y j + zli. We suppose that s is odd, to be justified under assertion[7l 

The first lattice point P of the set may be translated to the origin P = 1; now 
if fc = 1 the result becomes trivial. Otherwise since all distances squared ||P — Q\\ 
are integer, | US' — sP\\, so | \\S\\. By section [5] there exists X = GGD l{S, s), for 
which 

s = XX= \\X\\ =XX. 
Now the rational point S is transformed by X to the lattice: 

X = GCD l{S, s) = GCDz,(2s - 5, s) = GCDl(^, s) = GCDfi(S',s), 

so 

GCDi(S', s)^X and GCDfl(S', s) = X. 

Since \\S\\ > ||s||, the prime factorisation of S is no shorter than that of s; the extra 
factors S' must occur between X and X (for some factorisation), so 5 = XS'X. That 
is, S* is transformed to 

XSX = XXS'XX = s^S', 

where S' is integer. The scalar component of each side equals s^, so S' has unit scalar, 
and 5' e Z3. 

Also each lattice point Q is transformed by X to the lattice: 
Note that sQ — S — —{sQ — S) is pure, so 

s2 I ||sQ-S'|| =-{sQ~S)^ 

which is transformed to (scaled by s^) 

/ I \\XsQX -XSX\\ = -(XsQX -XSXf; 

so 

I sXQX-XSX, I XSX, \ sXQX 

via commutativity of scalars, and adding. That is, Q is transformed to XQX = sQ', 
where Q' is integer with unit scalar as before, and Q' ^1? . □ 

In contrast to dimension 2, permuting vertices may yield embeddings which are 
essentially distinct, in the sense defined in section [9l Consider the penultimate stage 
of the embedding procedure, where P,Q,R £ and only S remains to be located: 
since P,Q,R are essentially unique via assertion [51 at most 4 distinct configurations 
are possible, depending upon which vertex was initially specified as S. Therefore 
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Assertion 11. The number of essentially distinct embeddings of a Heronian tetrahe- 
dron into Tl? , as constructed via the GCD procedure, is at most 4- 

However, it will emerge in section ^ that there are embeddings which are not 
constructible in this fashion. Whether there is a bound on the total number of distinct 
embeddings is not known: modulo lattice isometry, it would have to be at least 36. 

Finally, complex numbers turn out to have been something of a waste of space all 
along: 

Assertion 12. Embedding in 1? is a special case of embedding in 1? . 

Proof. The claim is equivalent to the following. Given a rational point R = l + x/ri + 
2//r j in Q^, the quaternion GCD rotates it to some lattice point remaining in Z^; that 
is, X — GCDi(r' + a:;i + ?/j,r) = / + gk for some f,g£Z. To show this, 

GCDL(r + xi + yj, r) GCDi(r, {x mod r) i + (y mod r) j) 

remaindering modulo r; 

= GCDi(r, X mod r + [y mod r) k) 
via right product of second argument by — i (units are coprime); 

= GCD(7', X mod r + Ly mod r) 
setting k = t, its polar complex equivalent; 

= f + i^g, 

say. □ 

Notice how (polar) rotations cos t + l sin i e C correspond to cos t + sin t ■ k e H; 
whereas (axial) points x + ty G C correspond to 1 + xi + yj eH. 

9 Essentially Distinct Embeddings; Canonical Re- 
duction 

Here we consider only tetrahedra in K'^; generalisation to simplices in R" is a matter 
of routine, were it ever required. 

While there are infinitely many lattice embeddings of any Heronian tetrahedron, 
factoring out translations by fixing one vertex at the origin reduces this to a finite 
number. [For technical reasons however, it proves more convenient instead to further 
translate so that all coordinates are positive, with each component zero for at least 
one vertex.] 

Many of these fall into sets of (in general) 48 lattice isomorphs of one another: to 
eliminate these trivial variations from consideration, define the weak canonical form 
to be the earliest isomorph, modulo lattice symmetries, in numerical order — where 
coordinate components taken in the natural ordering are regarded as individual 'digits' 
in a 12-digit 'number'. 

However even with this refinement embeddings may proliferate, particularly (and 
paradoxically) for symmetric cases. They may be further reduced by a potential factor 
4 via permutation of vertices; so define the strong canonical form to be the earliest 
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weak form modulo vertex permutation. [The relationship between input edge lengths 
and vertices is disrupted, but can straightforwardly be recovered when required.] 

The symmetry types occurring in Heronian tetrahedra are of interest in this con- 
nection. The table records number of isomorphs for each group, along with the corre- 
sponding pattern of equal edges, and a (permuted) example where applicable. 
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The last two cases above are trivially impossible, since a unit equilateral triangle has 
irrational area •\/3/2; the 8-symmetric case is dispatched at [2] sect. 3 case 2(ii). 

The first Heronian triangle for which the lattice embedding is not also an axial 
pose is 

[30,29,5]; 

the first tetrahedron with some non- axial embedding, also the first without an axial 
embedding, is 

[160,153,25,120,56,39]; 
the first with some embedding not constructible via GCD from any axial pose is 

[888,875,533,533,875,888]; 

the first with 2,3,4 distinct embeddings via GCD are respectively 

[1073,975,448,495,952,840], 
[1360, 1092, 548, 975, 865, 663], 
[45100, 43911, 6929, 34476, 40544, 36975]. 

The particularly prolix case 

[8484, 6625, 6409, 6409, 6625, 8484] , 

yields 9(36) strongly(weakly) distinct embeddings, versus just 1(4) via GCD. 

10 Exhaustive Search for Lattice Embeddings 

The quaternion GCD-based construction establishes that some lattice embedding ex- 
ists for every Heronian tetrahedron case, and is furthermore fast — far more time is 
spent in reduction to canonical form than in construction. GCD embeddings have 
been computed for every case with diameter < 100, 000. [Indeed in 4 cases, new 
embeddings were found which had been overlooked — to considerable embarrassment 
on the part of its implementor — by the allegedly complete search described below.] 
We shall need the parameterisation generating all primitive solutions of the Dio- 
phantine equation a;^ + j/^ + = from the identity 

(p2 ^(^ -V? - v^ f + (2pu + 2qvf + {2vv - 2quf = (p^ + ^ v^f. (4) 
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[In [13], where it is credited to Mordell, this same idea is eniployed in an entirely 
different capacity — see end of section [T3l] 

Constraining the integer parameters p, q, u, v by 

i< P < [V^\ , 

[Vmax(0, u;/2 - p2)] < q < [y /w-p^J, 

< u < [■\/ min(p2 + q^jW — — q'^)\ , 
V = [■\/ w — p"^ — q^ — , 

before verifying +(7^ = w will ensure that a; < y < z; though some solutions 

may be imprimitive or duplicated. The complete solution requires merging scaled-up 
solutions for every divisor of w. 

Now the following algorithm can be designed to enumerate the set of all essentially 
distinct embeddings of a given case: 

1. Find some rational axial pose for vertices P,Q,R,S via section [7] noting that 
altitude of S has ambiguous sign; 

2. Generate complete solutions to x'^ + y'^ + z'^ ~ u'^,v'^ resp. via (|4]); 

3. With vertex P at origin, scan all possible lattice locations for Q, R, retaining 
partial embeddings where length QR^ = uP'] 

4. Solve over-determined homogeneous linear equations (via null space) for the 
quaternion X rotating face PQR from axial pose to lattice embedding; 

5. If X also rotates (either) S to the lattice, reduce the rotated pose P'Q'R'S' to 
canonical form and add to result set. 

A Magma implementation of this algorithm has currently found many more Z"^ 
embeddings for tetrahedron diameter < 20, 000, despite evidently remaining in need 
of more thorough development! A version employing assertion [12] was employed (after 
much ill-tempered adjustment of dimension) to verify uniqueness in 1? for the first 
1000 triangles. 

11 Worked Numerical Examples 

The following illustrate the complete procedure for computing a canonical lattice em- 
bedding of a Heronian tetrahedron via quaternion GCD applied to a rational axial pose. 
Throughout this section the quaternion s+pi + qj + rk will be denoted [s,p, g,r]. 
Input data comprise six edge-lengths. 

The first example involves the scalene case 

[2431, 2375, 1044, 2296, 2175, 1479]. 

Apply the currently selected permutation on vertices, say PQRS — > QRPS, and via 
section [7] construct the corresponding rational pose P,Q,R,S — 

[1,0,0,0], [1,1044,0,0], [29,18876,67925,0], [13,22620,8613,14616]; 

notice the odd scalar component (Cartesian LCD) in R, S. Since (inevitably) P, Q 
already lie in the lattice, via section |8] relocate R by setting 

X = GCDz,(i?, 29) = [-5, 0, 0, 2] 
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and transforming all four vertices hy X to P,Q^R,S = 

[1,0,0,0], [1,756,720,0], [1,-1144,2145,0], [13,10440,21837,14616]; 

then repeat to relocate S by 

X = GCDl(5,13) = [-2,2,2,1] 

transforming to P, Q,R, S — 

[1,0,0,0], [1,396,864,432], [1,396,-561,2332], [1,1740,783,1044], 

with every vertex now in the lattice. Finally translate to the positive octant; scan 
lattice symmetries and vertex permutations as per section [5J finding the canonical 
isomorph 

[1,0,0,396], [1,561,2332,0], [1,1344,1288,1740], [1,1425,1900,396]. 

The second example involves the same case, but permuting PQRS — > PQSR. Now 
the sequence follows rational pose 

[1,0,0,0], [1,2431,0,0], [13,17248,24360,0], [17,36575,13680,10260], 

hy X= [0,-3,-2,0] to 

[1,0,0,0], [1,935,2244,0], [1,2240,504,0], [17,26695,28500,-10260], 

by X = [0,3,-2,-2] to 

[1,0,0,0], [1,-1529,-1848,396], [1,-224,-1848,-1344], [1,-665,-2280,0], 

reducing canonically to 

[1,0,0,396], [1,224,1848,1740], [1,665,2280,396], [1,1529,1848,0]; 

different from the previous result and so representing an essentially distinct embedding. 
[In scalene cases there is no difference between strong and weak.] For this case all 24 
permutations of vertices yield one of these two results. 

Nonetheless, tedious exhaustive search (section (TU]) eventually discovers one more 
embedding, distinct from both above: 

[1,0,0,0], [1,224,672,2184], [1,665,1824,1368], [1,1529,1716,792], 

showing that not every embedding is related via GCD to an axial pose. 
The third example employs isohedral case 

[8484, 6625, 6409, 6409, 6625, 8484] . 

This embedding might be thought to lack variety: whatever the vertex permutation, 
any axial pose, say 

[1,0,0,0], [1,6625,0,0], [5,28224,31668,0], [5,4901,22932,21840] 

which happens to transform in a single shot hy X — [0, —1, —2, 0] to 

[1,0,0,0], [1,-3975,5300,0], [1,1680,8316,0], [1,3081,3536,-4368], 
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always yields the same strongly canonical 

[1,0,0,1401], [1,0,3016,7056], [1,0,8316,3081], [1,4368,4780,0]. 

Invoking only weak canonicity reduplicates this as 4 vertex-permuted variations. More 
noteworthy are the results of exhaustive search in this case, yielding in total 9(36) 
distinct embeddings. 

An intriguing feature of isohedral cases is the sparseness of their rotations: it 
appears that such X ^ \ have just two nonzero components, which furthermore are 
very small — eg. X = — (i-|-2j) above. For all 11 isohedral cases with diameter 
< 20, 000, there is some vertex permutation for which the rotations have height at 
most 4. 

12 Heronian Pentatopes in 

Consider the free pentatope PQRST specified by squared edge-lengths 

edge QP RP RQ SP SQ SR TP TQ TR TS 
square length 123232121 1. 

PQRST can be rationally posed in Q"*, with projective coordinates for its vertices: 

P= [1,0,0,0,0], 
Q= [1,1,0,0,0], 
i?= [1,0,1,1,0], 
^= [1,0,1,0,1], 
T= [3,0,1,2,2]. 

Taking the determinant shows PQRST is proper, with volume 1/24. 
However, 

Assertion 13. The simplex PQRST above has no lattice embedding in Z^. 

Proof. We may choose lattice axes so that T lies at the origin 

T= [1,0,0,0,0]. 

Since TS'^ = TR^ = TP'^ = 1 we may choose 

5= [1,0,0,0,1], i?= [1,0,0,1,0], P= [1,1,0,0,0]. 

Since QP^ = 1, and the volume is nonzero, the vacant column must be filled by 
choosing 

Q= [1,0,1,0,0]. 

Now RQ^ ~ 2, contradicting = 3 in the specification. □ 

The point of all this is that assertions Oil] require only that a polytope be rationally 
embedded, and have edge lengths squaring to integers, in order to be embeddable in I? 
or 1? [constraints garbled in the final paragraph of [13]]. The counter-example above 
demonstrates that these constraints are insufficient to ensure embeddability in Z^; 
it follows that any theorem concerning Heronian embeddability in higher dimensions 
must impose stronger conditions. 

That's always assuming that somebody can manage to find anything to embed. 
Sascha Kurz |10j recently completed the enumeration of primitive Heronian simplices 
for diameter < 600, 000. He reports 41563542 triangles, 2526 tetrahedra, and no 
pentatope at all — nought, zero, zilch. 
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13 Remarks on Alternative Approaches 



During the course of an interesting but in places unspecific paper covering many of 
the same topics as here, the proof of [8] Theorem 4.4 develops a construction for 1? 
triangle embeddings similar in spirit to our proof of assertioijH however, there are 
significant differences. 

In the first place no GCD is explicitly invoked: instead the LCD r is factored into 
rational primes p, then by induction each factor is dispatched via a new rotation X — 
obtained in essentially the same manner as above. From a computational viewpoint, 
this decomposition is both more elaborate and slower. 

When it comes to dealing with the set however, instead of each point dispatched 
in turn via a new rotation, every point is (partially) embedded via the same X . Com- 
putationally speaking, what is lost on the swings — the number of factors oi q — is 
(perhaps) regained on the roundabouts — the number fc of points. 

So we should enquire whether — in especially — our induction on k might also 
be aggregated, winning on both rides. Instead of transforming by X = GCDL(Q,g) 
for each point Q in turn, consider just setting X — GCDi(S', s) for the final point S 
of the rational pose, in the hope that X will also embed every other point. 

Clearly this could work only if all previous q\s. But the situation is actually 
rather worse: it turns out that, after rotation and reduction, the scalar Q' remains 
GCD(s/g, q) rather than unity. So rotated points may miss the lattice unless, for each 
scalar LCD q and prime factor p, either p divides q to the same power as p divides s, 
or p does not divide q at all. 

And in practice cases do frequently occur where, besides denominators involving 
various powers of the same prime, the rational pose has no LCD equalling their LCM: 
for example (for various permutations of the vertices) the early case 

[160,153,25,120,56,39] 

yields some axial poses with scalars g = 1, 1, 25, 65. Tough luck! 

Turning to R^, the paper provides a partial solution to the existence problem: [8] 
Theorem 4.6 asserts that computation has established the existence of suitable rotation 
matrices, provided all LCDs are products only of primes p < 37. Unfortunately no 
details are given: one conjectures that the method involves some variation of the 
'modularisation' employed elsewhere in attacking divisibility of volume, see section [2l 

A complete contrast to the complex approach to embedding in R2 is 

provided by 

|13j , where the problem is reformulated in terms of a quadratic Diophantine equation, 
for which the complete parametric solution Q is known. It seems improbable that an 
analogous attack could easily be mounted in R^: corresponding parameterisations for 
simultaneous quadratic and cubic Diophantine equations are (as far as we are aware) 
currently unknown, and quite possibly do not exist. 
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